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£_! ■ Abstract. 

o: 

We develop the fundamentals of a local regularity theory for prescribed Jacobian equa- 
'. tions which extend the corresponding results for optimal transportation equations. In this 

theory the cost function is extended to a generating function through dependence on an 
additional scalar variable. In particular we recover in this generality the local regularity 
theory for potentials of Ma,Trudinger and Wang, along with the subsequent development 
of the underlying convexity theory. 
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Oh! 1. Introduction 

< 

Let SI be a domain in Euclidean n-space, R™, and Y a mapping from x R x R™ into 
R™. The prescribed Jacobian equation, PJE, is a partial differential equation of the form, 

(1.1) det DY(-,u,Du) = ip(-,u,Du), 

where if> is a given scalar function on SI x R x R™ and Du denotes the gradient of the 
function u : SI — > R. Denoting points in SI x R x R n by (x, z,p), we see that the special 

case, 

(1.2) Y{x,z,p)=p, 
corresponds to the standard Monge- Ampere equation, 

(1.3) det D 2 u = ip(-,u,Du). 

We will always assume that the matrix Y p is invertible, that is det Y p ^ 0, whence we may 
write (1.1) as a general equation of Monge- Ampere type, 



(N 



^ (1.4) det[D 2 u- A(-,u,Du)] = B(-,u,Du), 

where 

(1.5) A = -Y-\Y X + Y z ®p), B = (det F p )~V- 

A function u G C 2 (fl) is degenerate elliptic, (elliptic), for equation (1.4), henceforth called 
admissible, whenever 

(1.6) D 2 u- A(-,u,Du) > 0, (>0), 

in Q. If u is an elliptic solution of (1.4), then the function B(-,u, Du) is positive. Ac- 
cordingly we will assume throughout that B is at least non-negative in Q x R x R™, that 
is if) and det Y p have the same sign. 

The second boundary value problem for the prescribed Jacobian equation is to prescribe 
the image, 

(1.7) Tu(p) := Y(-, u, Du)(n) = SI*, 
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where Q* is another given domain in M. n . When ip is separable, in the sense that 
(1.8) \ip(x,z,p)\ = f(x)/goY(x,z,p), 

for positive /, g G L 1 ^), L 1 ^*) respectively, then a necessary condition for the existence 
of an admissible solution, for which the mapping T is a diffeomorphism, to the second 
boundary value problem (1.1), (1.7) is the mass balance condition, 

(1-9) f f = f 9- 

For the standard Monge- Ampere equation (1.3), with Tu = Du, the classical solvability of 

the second boundary value problem, under the mass balance condition (1.9), was proved 
by Delanoe, (n = 2), [4], Caffarelli [2] and Urbas[29], under the hypothesis that both 
domains, Q and Q* are uniformly convex. As already pointed out in [12], (1.7) implies a 
nonlinear oblique boundary condition. A weaker interpretation of the boundary condition 
(1.7) arises through optimal transportation, in which case Caffarelli [1] proved that the 
convexity of the target Q* suffices for local smoothness of solutions. 

Interest in the general case was stimulated in the last decade through its application 
to regularity in optimal transportation, ([19, 26]). Here we are given a cost function 
c:R"xlMK and the vector field Y is generated by the equation, 

(1.10) c x {x,Y{x,p))=p, 

which we assume to be uniquely solvable for p G c x (Q x Q*), with non- vanishing determi- 
nant, that is det c XjV (x, y) ^ 0, for all (x, y) G f2xf2*. In the corresponding Monge-Ampere 
equation (1.4), we have 

(1.11) 

A(x,z,p) = c xx (x,Y(x,p)), B(x,z,p) = det c XtV (x, Y(x,p)tp(x, z,p). 
Conditions for local regularity were found in [19] and for global regularity in [26], with 

the latter being extended to general prescribed Jacobian equations in [21, 22]. 

In this article we consider a more general situation where the cost function is replaced 
by a smooth generating function G : M. n x M. n x K. — > R, with the resultant vector 
field Y also depending on u. This enables the corresponding local theory to embrace 
recent work in near field optics, [10, 11]. The assumptions on the generating function G 
are parallel to those introduced for cost functions in optimal transportation in [19, 26]. 
Accordingly we let U be an open set in W l x W 1 and I an open interval in R. For points 
(x,y) G U, we denote their corresponding projections by U* = {y G M, n \ (x,y) G U}, 
U y = {x G R n I (x, y) G U} and write = \j{U y \ y G M n } and = {j{U* \ x eW 1 } 
Denoting points in I by z, we assume that G is smooth in U x 7, G z ^ and 

Gl: For each (x, y) G U, there exists an open interval I(x, y) C I such that the mapping 
(G x , G)(x, •, •) is one-to -one in y G U*, z G 7(x, y), for each x G US 1 ^ . 

G2: For each (x,y) EU,z G I(x,y), det E(x,y, z) ^ 0, where E is the n x n matrix 
given by 

(1.12) E = [E Xty \ = [G XtV - {G z )~ l G XtZ ® Gy]. 

From Gl and G2, the vector field Y, together with a scalar function Z, are generated by 
G through the equations, 
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(1.13) G x (x,Y,Z)=p, G(x,Y,Z)=u. 

The significance of the additional function Z will become apparent later. Note that the 

Jacobian determinant of the mapping (y, z) — > (G x , G)(x, y, z) is G z det E, 7^ by G2, so 
that Fand Z are accordingly smooth. Also by differentiating (1.13), with respect to p, 
we obtain Y p = E~ x . Next using (1.5) or differentiating (1.13) for p = Du, with respect 
to x, we obtain that the prescribed Jacobian equation (1.1) is a Monge- Ampere equation 
of the form (1.4) with 

(1.14) A(x,u,p) = G xx [x,Y{x,u,p),Z{x,u,p)\, 

B(x, u, p) = det E(x, Y, Z)ip(x, u, p) 

and is well defined in domains Q G US^> for Y eU*,Z e I(x,Y),x G Q. Note that the 
latter restrictions may automatically place constraints on u and Du. 
In the optimal transportation case, we have 

(1.15) G(x,y,z) = c(x,y) - z, G z = -1, I = I(x,y) =R, 

and we recover (1.10) and (1.11) above. 

Note also that by adjusting the dependence of G on z, we can fix the interval /, as well 
as the sign of G z , as we wish. For convenience with other examples, we will assume either 
/ = (0, 00) or I = R as above and assume G z < 0. Let us also denote 

r = T(U) = {(x,y,z) eUxI\ze I(x,y)} 

We mention also that the idea behind conditions Gl and G2 is to determine the mapping 
T from a tangential intersection of the graphs of the functions u and G(-, y, z). If the graph 
of G is also a local support from below, we obtain the ellipticity condition (1.6), which 
corresponds to a local convexity. The geometric picture is further amplified in [23] , where 
the theory developed here is already outlined. Knowledge of the defining function G can 
also lead to a more efficient derivation of the Monge- Ampere equation, using (1.14) rather 
than computing it directly from (1.5). Also note that when the graph of G is a local 
support from above we obtain the complementary ellipticity condition, 

(1.16) D 2 u — A(-, u, Du) < 0, 

corresponding to a local concavity. In this case by setting u~ = —u,A~(-,u,p) = 
—A(-,—u,—p), we also obtain a degenerate elliptic equation of the form (1.4) for the 
function u~. Furthermore if n = 2 and B > in (1.4), then either of the strict versions 
of (1.6) or (1.16) must hold. 

For the regularity of weak solutions we need a global convexity theory. To develop this, 
we in turn need a dual condition, which in the optimal transportation case is simply 
obtained by interchanging x and y. Using the property that the generating function G is 
strictly monotone with respect to z, we introduce a dual function G* = H on U x R by 

(1.17) G[x,y,H(x,y,u)] = u. 

Clearly H is well defined whenever (x,y) G U and u G G(x,y, •)(/). Furthermore we 
have the relations 
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(1.18) H x = —G x /G z , H y = —G y /G z , H u = 1/G Z . 
This motivates the following dual condition: 

Gl*: The mapping Q := —G y /G z is one-to-one in x, for all y G U^ 2 \z G I(x,y). 

Writing J(x,y) = G(x,y, -)I(x,y), we see that condition Gl* corresponds to Gl with x 
and y interchanged and I(x,y) replaced by J(x,y). Furthermore the Jacobian matrix of 
the mapping x — > Q(x,y,z) is —E/G z so its determinant is automatically non-zero when 
condition G2 holds. Analogously to T above, we may also denote the dual set, 

T* = {(x,y,u) G U x R | u G J(x,y)}. 

Our next conditions extend the conditions A3 and A3w introduced for regularity in 
[19, 26] and are expressed in terms of the matrix function A in (1.4), namely: 

G3 (G3w) A^VkVi ■•= (D PkPl A l3 )UjVkVi > (>) 0, 

for all (x, Y) eU,Z e I(x, Y), f , rj G R n such that = 0- 

As in [21], we may write equivalently that A is strictly regular, (regular), in the set 

{(x,u,p) | x G U {1) ,u = G(x,y,z),p = G x (x,y,z),y G U*,z G I(x,y)}. 
For our convexity theory we will also need a condition that the matrix function A is 

monotone with respect to u, that is 

G4 (G4w) /Vt, ; C<, > (>) 0, 

for all (x, Y)eU,Z G I(x, Y),£ G W l . 

We illustrate the above conditions with the example of a parallel beam from [17]. For 
U = R n x R n and I = (0, oo), we define 

(1.19) G[x ^^ z) = L_ z _\ x _y\\ 
Then G satisfies G1,G2,G1*,G3, G4 with 

G z = -^{z~ 2 + \x-y\ 2 ) <0, G x = z{x-y), G y = z{y-x), 

E = z{(l + z 2 \x - y\ 2 )I - z 2 (x -y)®(x- y)}/(l + z 2 \x - y\ 2 )), 
detE = z n (l - z 2 \x - y\ 2 )/{l + z 2 \x - y\ 2 ) >0,z e I(x,y), 

1 1 

^ = — i 7 2 i i f2u72> y) ~ (0) i r),^(^y) = (o,oo). 

■u + (u 2 + \x — y\ l Y IA \x — y\ 



1 - \Du 



The corresponding Monge-Ampere equation, 

(1.20) det(i) 2 M+ (1 ~i D " |2) j}= r 

v ; 1 2m J (1 + |Dm| 2 )(2m)" 

is well defined for w > and \Du\ < 1. 



2^n+l 
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In the next section, we show that a convexity theory, with respect to generating func- 
tions, replicating usual convexity, and more generally the optimal transportation case as 
developed in [9, 13, 25, 27, 32] can be built under conditions Gl, G2, Gl*,G3w and G4w. 
In particular we prove that the local convexity of smooth functions implies their global 
convexity for appropriately convex domains, that normal mappings are determined by 
sub-differentials and that sections and contact sets are convex in the extended sense. As 
well we prove that domain convexity with respect to generating functions is a special case 
of that determined by vector fields in [21, 22]. Our proofs follow the optimal transporta- 
tion presented for example in [25] , although they also depend on more intricate 
calculations. 

In Section 3, we derive the dual Monge- Ampere equation, satisfied by the G-transform 
v of a G-convex solution u of (1.1), given by 

(1.21) v{y) = u* G (y) = sup n H(-, y, u). 

and prove that conditions G3 and G3w are invariant under duality. 

Section 4 is devoted to the existence and regularity of generalized solutions of the 
second boundary value problem (1.7), (1.8), with initial and target domains satisfying 
f2 x f2* c U. For existence we follow the approach in [3, 10], which corresponds to the 
existence of potential functions in optimal transportation. Here we need an additional 
condition to control gradients of solutions, which we may express as: 

G5: There exists constants m > —oo,K Q > 0, such that (m ,oo) C J(x,y) and 

\G x (x,y,z)\ < K 
for all x G f2, y G G(x, y, z) > m 

The example (1.18) clearly satisfies G5 for m = 0,K = 1, while in the optimal trans- 
portation case (1.15), m = — oo. Note that by translation, we may always assume either 
mo = or m = — oo. Under condition G5, there exists a generalized solution u with 
u(xq) = Mo, provided Uq > m + K d(xo,dQ). More generally we can also replace G by 
composites fJ>(G) for suitable smooth functions ji G C 1 (m , oo), // ^ 0. This enables us 
to embrace more examples such as near field reflection from a point source, as treated in 
[10, 11]; (see Section 5). 

Following [19], we then prove local regularity under conditions G1,G2,G1*G3 and G4w 
for target domains satisfying the appropriate convexity conditions. We remark that the 
existence of globally smooth elliptic solutions under corresponding conditions, including 
stronger domain convexity conditions, follows from the theory of the general prescribed 
Jacobian equation in [22] and its extension to near field reflection problems in [17]. 



2. Convexity theory 

We begin with the appropriate definitions of convexity with respect to generating func- 
tions. 

Let O be a bounded domain in U^ l \ denote = (~] x( znU* and let G be a generating 
function on U x /, satisfying conditions Gl and G2 , with / an open interval in R as in 
the previous section. A function u G C°(fl) is called G- convex in Q, if for each x G f2, 
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there exists y G and z G I(x ,y ) such that 
(2.1) u(x Q ) = G(x Q , y Q ,z Q ), 

u(x) > G(x,y ,z ) 

for all isSl. If u is differentiable at xo, then y = Tu(xo) = Y(xo, u(xo), Du(xq)), while 
if u is twice differentiable at x , then 

D 2 m(x ) > DlG(x ,y ,Zo), 

that is m is admissible for equation (1.4) at x . If u G C 2 (Q), we call u locally G- convex 
in f2 if this inequality holds for all Xq G f2. We will also refer to functions of the form 
G(-, yo, Zq) as G-affine and as a G-support at x if (2.1) is satisfied. 

As in the optimal transportation case [19], we also have corresponding definitions of 
domain convexity. There are various possibilities depending on what quantities are fixed. 
For our purposes here we make the following definitions: 

the domain Q is G- convex with respect to y G Wq, Zq G / if the image Qo(Q) := 
—Gy/G z (-,y ,z )(£l) is convex in M n ; 

the domain Q* is G*-convex with respect to (xo,Uo), where Xo G f2 and u G J(xo,y) for 
all y G Q*, if the image Po(^*) : = G x [x , -,H(x , -,u )](fi*) is convex in IR". 

Note when we use condition G3w for convexity results and their consequences, we will 
assume at least that the convex hulls of the images Qo(£l) and P (^*)) he in Q(T), G X (T) 
respectively. 

The second definition is clearly a special case of the notion of y*-convexity in [21], as 
it is equivalent to the set Pq(Q*) = {p G M ra | Y(xq,Uq,p) G being convex. Moreover 
we will define the domain Q* to be G*-convex with respect to a function u G C°(Q) if Vt* 
is G*-convex with respect to each point on the graph of u, that is the sets P X (Q*) = {p G 
M n | Y(x,u(x),p) G il*} are convex for each x G f2 The relationship of the first definition 
and the notion of F-convexity in [21] is treated in Lemma 2.4 below. Our main Lemmas 
2.1,2.2 and 2.3 extend corresponding results for optimal transportation in [27, 9, 14] and 
[5]. 

Lemma 2.1. Assume Gl,G2,Gl*,G3w and G^w hold in 11 x I and that u G C 2 (Q) is lo- 
cally G- convex in Q. Then if Q is G- convex with respect to each point in (Y, Z)(-,u, Du) (Q), 
u is G-convex in Q. 

More generally, we have for any x G Q,yo,z = Y, Z(xq, u{xq), Du(xo)), the G-affine 
function G(-,yo,z ) is a G-support, provided Q is G-convex with respect to (yo, ^o)- 

We defer the proof of Lemma 2.1, until the end of the section and proceed now to 
consider the corresponding extensions of normal mappings and sections. 

Let u G G°(f2) be G-convex in Q. We define the G -normal mapping of u at x G fl to 
be the set: 

Tu(x ) = {y eUn \ u(x) > G(x,y , H(x Q ,y ,u f) for all x G fij, 

where Uq = u(xo). Clearly Tu agrees with our previous terminology when u is differen- 
tiable and moreover in general 



Tu(x ) C Y(x ,u(x ),du(x )), 
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where du denotes the sub differential of u. 

Lemma 2.2. Assume Gl,G2,Gl*,G3w,G4w hold in U x I and suppose u G C°(Q) is 
G-convex in Vt. Then for any x G fl, we have 

Tu(x ) = Y(x ,u(x ),du(x )). 

For a fixed yo G Tu(xq), corresponding zo = H(xo,yo,uo) and a > 0, we define the 
G -section, S a by 

5 CT = S a (x ,y ) = \x G Q | < y ,^o - c)} 
and lower contact set So by 

S = S (x ,y ) = {i G O u(a;) = G(x,y , z )]. 

Lemma 2.3. Assume Gl,G2,Gl*,G3w,G4w hold in U x 7 and suppose u G C°(f2) «s 
G-convex in Q, with Q itself being G-convex with respect to yo G Tu(xq) and zq = 
H(xo,yo,u(xo)). Then the sets S a and So are also G-convex with respect to (yo,z ). 

The proofs of Lemmas 2.1,2.2 and 2.3 reduce to calculations that are more complicated 
versions of those which underwrite our starting point in [25, 26], relating domain convexity 
to the Monge- Ampere equation (1.4) . We will adopt the following notation: 

d d d d 

(2 ' 2) Gi > " > z >- ■ dx~'"dy~ "'ih"'7h (: - 



E — E 



E 



= E- 1 = 



so that in particular we have the following formulae for differentiation with respect to the 
p variables: 

G 

(2.3) Z p = —Y p — 



P G Z 

D p = Y p (D y - ^D z ) 

n = E r ' k (D v — — D. 

Pk y yr Q ' 

n p — p r ' k ( p - p 



G._ 

Next from condition Gl*, we infer the existence of a mapping X, defined uniquely by 

G„ 



(2.4) -g-(x(y,z,q),y,^ = -q 

G 

for all q G —-gr(-,y, z)(Q). It follows that 

(2.5) X q = -G Z Y P , E = -G Z X- 1 
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while for fixed y, z, we have the following formulae for differentiation with respect to the 
q variables: 

(2.6) D q . = —G z E i ' r D Xr 

D 2 qm = C: I^CCjI),,, + G z E^D Xs (G z E^)i^D Xr 

= G 2 z E hr E j ' s ^jl D XrXs + — E k ' l (G z G rs ^ k — G rS)Z G^D Xl 



G z 

= G 2 z E l ' r E j,s £ i £ j ^D XrXa - (p pi G r ^D xl - — G s>z D Xr ^ 

= G 2 z E^E^\D XrXa - (D pi G rs )D x ^ + ■lE-X,G,,n, lL . 

for any unit vector, £ G W 1 and q^ = q.£. 

The proofs of Lemmas 2.1, 2.2 and 2.3 follow from formulae (2.6). First we note also 
that our definition of G-convex domain aligns with that determined by the vector field Y 
in [21]. 

Lemma 2.4. Assume G1,G2 and Gl* hold inhi x I with dVt G C 2 . Then Q is G -convex 
with respect to y <EU*,z G I(x, yo) for all x G Q if and only if 

(2.7) [Di 7j(x) - G ij:Pk (x, y , z ) ^ k (x)} n tj > 0, 
for all x G dQ, unit outer normal 7 and unit tangent vector r. 

We prove Lemma 2.4 by applying formula (2.6) to the distance function d = dist(-, dtt) 
and using the orthogonality of Dd = —7 and r to cancel the last term in (2.6). 

To prove Lemmas 2.1, 2.2 and 2.3, we let u be locally G-convex in Q and for a fixed 
point x ,y <EU,z = H(x ,yo,u(x )), define the height function, 

(2.8) h(x) = u(x) - G(x,y ,z ). 

Making the coordinate transformation, x — > q — — G y /G z (x, yo, Zq), and setting y = 
Tu{x) = Y(x,u(x), Du(x)), z = Z(x,u(x), Du(x)) we compute, using formulae (2.6), 

(2.9) 

D q( h = -G z E^ k ^D Xk h, 

D 2 (q( h = G 2 z E l > k E^[D XkXl - D Pr G kl D Xr }h&j + 2/v'Y/,.,/V'^, 

> G 2 z {G m (x, y, z) - G ki (x, y , z ) - D Pr G k i(x, y , z ) 

[G r (x, y, z) - G r (x, yo, Zo^E^E^j + 2E"G 9 , l 4 j D qe h, 

for any unit vector, £ G R n . Setting u = G(x,y ,z ),p = G x (x,y,z),p = G x (x,y ,z ) 
and using condition G4w, we then have for h < 0, 

(2.10) 

Dg^h > G 2 z [A k i(x,uo,p) - A k i(x,uo,Po) - D Pr A k i(x,u ,p ) 
(p r - porW^E^ + E^G S;Z CjD q( h, 

> l -D PrPs A kl {x,uo,p*)D r hD s h^)E l ' k E^ 1 ^ + E^G^D^h, 

> -K\D qe h\, 
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by condition G3w and Taylor's formula, for some p* on the straight line segment £ joining 
p and po and constant K depending on G,Q and £. Setting q = q(x ,y , z ), it follows 
then that if D v h(x ) > where rjj = E h:i (q i — q 0i ), we have the inequality, 

(2.11) < D v h(x ) < Ch(x) 

for some positive constant C, depending on G,Q and Tu(Q). This result thus extends the 
fundamental lemma in [25]. By taking yo = Tu(xq), we conclude Lemma 2.1. 

To prove Lemma 2.2, we first note that a G-convex function u is semi-convex so that 
at any singular point Xo, its subgradient du(xo) is a convex set whose extreme points are 
limits of points of differentiability. The result then follows by showing that the image 
-Po : = G x [(xo, -,H(xo, •, uo)]Tu(xo), where uo = u(x ), is convex in W 1 , (that is Tw(xo) is 
G*-convex with respect to (xo,wo))- Accordingly we fix two points 2/1,2/2 G Tu(xq) and 
define corresponding G-afnne functions, Ui(x) = G(x,yi,Zi) for Zi = G(x ,yi,u ),i = 1,2. 
Then for any p , lying on the straight line segment joining Dui(x ) and Dm 2 (jo), and 
y = Y(x ,u ,po) , either u\ or u 2 satisfies the above condition, D^h^Xo) > 0, with 
respect to a fixed x G Q. Consequently G(x, y , z ) < u(x) and we conclude y e Tu(x ) 
as required. 

Lemma 2.3 follows immediately from the differential inequality (2.10), if u e C 2 (Vt). 
Otherwise, we may argue along the lines of [25, 28]. Assuming S a CC for any point 
xi G dS a ,ui = u(x\),yi G Tu(xi),zi = H(xi,yi,Ui), the inequality, 

(2.12) G(x, yi,z ± ) < G(x, y ,z - a), 

holds for all x G S a . Making the transformation x — > q = Q(x,yo,zo) as above, we need 
to show the transformed domain Qq = Q(-, yo, Zo)(Q) is convex. If Qq is not convex, there 
must be a straight line segment i joining two points in Q , containing a boundary point 
qi = q(xi). Applying the differential inequality (2.10) along £, we reach a contradiction 
with (2.11). 

We remark that Lemma 2.3 is equivalent to sub-level sets of the function —G*(-,yo,u) 
being convex with respect to the transformed coordinates q. In the optimal transportation 
case, —G*(-,yo,u) = u — c(-,yo), and this result is proved and used in [5, 14]. 



3. Duality. 

In this section we derive the dual Monge-Ampere equation for the G-transform v of 
an elliptic solution u of the prescribed Jacobian equation associated with a generating 
function G and prove the invariance of conditions G3, G3w under duality. Let u G 
C 2 (Q) be elliptic for (1.4) (1.14) and suppose that the mapping Tu defined by (1.7) is 
a diffeomorphism from to Q, to the target domain where G satisfies G1,G2,G1* and 
Vt x Vt* C U. Then we define the local G-transform of u by 

(3.1) v(y) = u* Gjloc (y) = H(T^y,y,uoT-\y)) 

where H = G* is the dual generating function introduced in (1.17). When u is G-convex 
in f2, then u* Gloc agrees with the G-transform defined in (1.20) . From (3.1) we have, for 

yen*, 
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(3.2) 

and hence 



Dv(v) = -% l {t- 1 v,vAv)) 



G 



(3.3) 



T- 1 = X(-,v,Dv) 



by Gl*, where X is defined by (2.5). Consequently if u is an elliptic solution of the pre- 
scribed Jacobian equation (1.1), (1.19), then v is an elliptic solution of the dual prescribed 
Jacobian equation, 



(3.4) 
where 



det 



D 2 v — A*(-, v, Dv) = B*(-, v, Dv) 



(3.5) 



A*(y,z,q) = H yy [X,y,u(X)] 



^2/ y 



G, 



{X ) y ) z)®q 



B*(y,z,q) = 



detH, 
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x,y 
1 ^ n 



G; 



foX 
detE 



foX' 



satisfying the dual second boundary value problem: 



(3.6) T*v(Q*) :=X(;v,Dv)(£l*) = Q. 

We can then formulate the dual conditions, 
G3*(G3*w) DwAUZ&rwt >,(>)0, 

for all (X, y)eU,ze I(X, y), £, 77 e M n , £ • V = • 

Theorem 3.1. Conditions G3, (G3w) and G3* , (G3*w) are equivalent. 

Proof : It will be convenient where possible to express our formulae in terms of the 
matrix function [E it j] defined in Condition Gl and the vector function [Qj] given by 

Q,i = Qi = ~ m Condition Gl*. Partial derivatives of these quantities will then align 
with our notation in (2.2). As already indicated in Section 1, we have Qij = —Eij/G z . 
Accordingly we can express the last formula in (2.3) in the form 



(3.7) D Pk Aij — E ' (^Eij r Qj,rGi^ z ^J 

= E r > k (E ijir + -J^GiJ) 

= E r,k E ijiT . + -if- 5 i>k . 
Now differentiating, with respect to pe, we obtain 
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(3.8) 



D 2 pkPe A t3 = (D Pe E r > k )E t3 , r + E^ k D Pe {E lhr ) + D n &±)8 jk 



G, 



so that using (2.3) again we have 



(3.9) A k j — E r ' k E s ' e ^Eij )rs + Eij y r yZ Q yS — E r ' s Eijy[E s ^ rs + E s >^ z Q iS ^J j 

G 7 



To get the corresponding formula for A*, we first write, 

1 / E r fa \ 

Qr,k£ = -Tr[E r ,ke — pT-C^J 



and 



Then we have 



^Cr,k,z \ JZ/ r,k,z ^ zz J 



(3.10) D qt A* kf = E l ' r [E r ^i + E rtk:Z Q t ij — — (p,t,z + G zz Q^5 ik 

+ Q,k,z8i£- 

Consequently we obtain the dual formula, 



(3.11) 



D 2 qiqj A* ki — —G z E t,r E j,s ^E rStki + Eij tk>z Q t £ — £ r ,s E rs>r >(E s > >k £ + E s > >ktZ Q t ij j 



+ E t,r Er tkjZ 5j£ — D, 



<h 



G 7 



G : £ jZ + G ZZ Q,£) 



5ik + D qj [Q,k, 



Formulae (3.9) and (3.11) will hold for all (x,y,z) G T and imply the equivalence of 
conditions G3(G3w) and G3*(G3*w), using Condition G2, (that is det£ ^ 0), and the 
orthogonality of £ and rj. 



4. Existence and regularity. 

We begin with the definition of generalized solution. Let Q and Q* be bounded domains 
in R n , with Q x Q* cU, and let u G C°(fi) be G-convex in fi, with conditions G1,G2, 
Gl* satisfied. Following [19], we introduce a measure \x = fi g [u] on fi, for g > G L 1 ^), 
such that for any Borel set E C fi, 



(4.1) /i(E) = jf 



5 

Tu(E) 

To prove that // is a Radon measure, we can extend the argument in [19], using the 
G-transform v, defined in (1.21). Note that v will be G*- convex in Q*. Defining the dual 
G*— transform on Q*, v*, by 
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(4.2) v*(x) = sup G(x,y,v(y)), 

yen* 

we obtain v * = u in f2 and furthermore y G Tu(x) if and only if x G T*v(y), where T* 
denotes the G*-normal mapping on Q*. Since the functions u and v are semi-convex and 
hence twice differentiable almost everywhere in Q and Q*, respectively, we then infer that 
the subset of points y G Q*, such that y G Tu{xi) fl Tu(x 2 ) for some X\ ^ x 2 , G f2, 
has measure zero and from this it follows that ji is countably additive. Also extending 
the argument in [19], we have that ji is weakly continuous with respect to local uniform 
convergence, namely if {u m } is a sequence of G-convex functions in Q, converging to u, 
then the sequence of measures {/i s [w m ]} converges to fi g [u] weakly. A G-convex function 
u on Q is now defined to be a generalized solution of the second boundary value problem 
(1.7) for equation (1.4), (1.14), under the mass balance condition (1.9), if 

(4.3) fi g [u] = Vf 

where Vf = fdx and g is extended to vanish outside Q*. More generally we can replace 
Vf by any finite Borel measure v on Q. We remark that this notion corresponds to that of 
generalized solution of Type A in [11], where a generalized solution of Type B corresponds 
to the dual notion, 

(4.4) //>](£*):= / lt f = S(E) 

for any Borel set E* C that is fJ,/[v] = u*, where v* = gdy. As in the optimal 
transportation case in [19], the two notions are equivalent provided the measures v and 
v* have densities / G L 1 (fi), g G L x (f2*) respectively. Therefore we have: 

Lemma 4.1. Let u be a generalized solution of the second boundary value problem (1.7) 
for equation (1.4), (1.14), under the mass balance condition (1.9), and let v be the G- 
transform of u. Then v is a generalized solution of the second boundary value problem 

(3.4) , (3.6). 

Next we formulate an existence theorem under condition G5 which extends the optimal 
transportation case. 

Theorem 4.2. Let il and f2* be bounded domains in M n , with Q x Q* C U and let G 

be a generating function satisfying G1,G2,G1* and G5. Suppose that f and g are positive 
densities in L 1 (f2) and L l {Vt*) satisfying the mass balance condition (1.9). Then for any 
xq G fl and Uq > m + K\, where Ki = i^ dist(x , dQ), there exists a generalized solution 
of (1.4), (1.7) satisfying u(xo) = u . 

The proof of Theorem 4.2 follows the approach in [3, 10] using approximations by 
piecewise G-affine functions to solve the dual problem (4.4), where v* is approximated by 
a linear combination of Dirac delta measures. Condition G5 implies that a generalized 
solution u whose graph passes through (xo,Uq) satisfies u > mi for some mi > m , 
together with an a priori bound: 

(4.5) |«|o,i < K, 

with constant K depending on K Q ,uo and diam Q. The solvability of the dual problem, 
with an arbitrary finite measure v* then follows from the weak continuity of fi*. Theorem 
4.2 then follows from Lemma 4.1. 
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Taking account of our remark following the formulation of G5 in Section 1 and choosing 
fj, = log with m = 0, we see that if we replace the gradient bound in G5 by 

(4.6) \G x (x,y,z)\<K G(x,y,z), 

we conclude that there exists a generalized solution u such that u{xq) = u for any given 
Uo > 0. We remark also that instead of prescribing u(xo) in Theorem 4.3, we may fix 
instead u min — m 1 > m . 

The theory developed in the preceding sections can now be applied to prove interior 
regularity results for generalized solutions. A crucial lemma for this purpose is Lemma 2.2 
on the characterization of the mapping Tu by the sub differential du, which shows that 
the concept of generalized solution is local, that is it is the same for subdomains.This 
property was originally overlooked in [19] but rectified subsequently in [27], with a more 
direct approach found in [9]. In this paper we will just prove interior smoothness under 
condition G3 following the approach in [19]. Further extensions of optimal transportation 
regularity depend also on Lemma 2.3 on the G-convexity of sections but we will not pursue 
these here. First we state some more lemmas which extend the corresponding results in 
[19]. 

Lemma 4.3. Let u be a generalized solution of (1.4), (1.14), (1.7). Suppose f > in Q 
and fl* is G-convex with respect to u. Then Tu{Q) C Vt . 

Lemma 4.4. Let Vt be a bounded domain in W 1 and u, v e C°(Q), G-convex functions 
satisfying u > v in Q and u = v on dQ. Then Tu(Q) C Tv(Q). 

The proofs of Lemmas 4.3 and 4.4 are essentially identical with the optimal transporta- 
tion case in [19]. 

Now we formulate the interior regularity result, which extends the main result in [19]. 

Theorem 4.5. Let u be a generalized solution of (1.7) with positive densities f G 
C hl (Q),g e C 1A (Q*) with /, 1// e L°°(Q),g, 1/g e L°°(Q*) and with generating function 
G satisfying conditions, G1,G2,G1*,G3 and G^w. Suppose that Vt* is G* -convex with re- 
spect to u. Then u e C 3 (f2) and is an elliptic solution of (1.4), (1.14). Furthermore if Q 
is G-convex with respect to v = u* G , then Tu is also a diffeomorphism from Q to Vt*, with 
v an elliptic solution of the dual equation (3.4). 

To prove Theorem 4.5 by the method in [19], we need to solve approximating Dirichlet 
problems, 

(4.7) det[D 2 w - A(-,w,Dw)] = B(-,w,Dw) in B r , 

w = u m on dB r , 

where B r is a small ball in Q and {u m } is a sequence of smooth functions converging 
uniformly to u. Specifically, we set 

U m — Uh m ~ C\x\ 2 , 

where u = u + C\x\ 2 is convex in B 2r and Uh m denotes the mollification of u with h m — > 0. 
The essential difference here with the corresponding equations in [19] is the dependence 
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of the matrix function A and scalar function B on u. However for for sufficiently large k 
and small r, the functions 

v = v m = u m + k(\x\ 2 - r 2 ) 

will still be sub barriers for (4.7) and again for sufficiently small r, we also have a classical 
comparison principle from which solution and gradient estimates can be inferred as in [19] . 
Global and interior second derivative estimates now follow as in [19] taking account of 
Remark 4.1 there; (see also [21, 22, 26]) and we conclude the unique classical solvability of 
the Dirichlet problem (4.7), together with that of the limiting problem for w G C°' 1 (B r ) n 
C 3 (B r ) satisfying w = u on dB r . Now using the monotonicity Lemma 4.3 and Lemma 
2.2, we infer w = u in B r and hence the smoothness of u in Q. 

We point out however that we have implicitly used a stronger condition on the target 
density g in the above argument, namely g G G 1,:l (Jl ). To use only the local smoothness 
of g we need to prove in advance that u G C 1 (fi), (to keep T(B r ) away from dQ*), or 
equivalently by duality that u is strictly G-convex. This may be accomplished under the 
hypotheses /, 1/f G L°°(Q),g,l/g G L°°(fi*) as in [27] or [13]. Moreover, as in [13], we 
also infer u G G 1,a (f2) for some a > 0. 

From the Schauder theory, we also obtain that if / G C°°(Q),g G C°°(Q*), then u G 
G°°(f2). As in the case with optimal transportation regularity results may be refined 
when perturbation arguments using Lemma 2.3 are employed instead of Lemma 2.2. 
In particular, when we replace C 1 ' 1 by C°' a in Theorem 4.5, we obtain from [18] that 
u G C 2 (Q). It would be interesting though to have extensions to G3w along the lines of 
[5, 30]. In this respect we note that Lemma 2.3 does not imply the G-convexity of level 
sets of a normalised potential u in the general case. 

We also remark here that if condition G3w is violated then Lemma 2.2 is no longer true 
in general and that generalized solutions are not necessarily continuously differentiable. 
This follows in the same way as the optimal transportation case in [13, 19]. Similarly from 
[19], the G*- convexity of the target Q* is also necessary for regularity. An interesting 
consequence is then that if the generating function satisfies G3, then generalized solutions 
of the complementary problem associated with condition (1.16) cannot be continuously 
differentiable in general. 

Combining Theorems 4.2 and 4.5 we get an existence result for locally smooth solutions. 

Corollary 4.6. Let f2 and il* be bounded domains in M. n , with f2 x f2* dlA and let G be a 

generating function satisfying Gl,G2,Gl*,G3,G4w and G5. Suppose that f G C 1 ' 1 ^ G 
C 1 ' 1 ^*) with f, 1/f G L°°(Q),g,l/g G L°°(Q*) and that f and g satisfy the mass balance 
condition (1.9). Then for any x G Q and u > m + K\, there exists a G-convex, elliptic 
solution u G G 3 (f2)o/ (1.4), (1.7) satisfying u(x ) = u and Tu(Q) = Q* a.e., provided fl* 
is G* -convex with respect to Q x (mo, oo). If also Q is G-convex with respect to Q* x /, 
then Tu is a diffeomorphism from Q onto Q* . 

We remark that the interval (mo-oo) in the G*-convexity hypothesis may be replaced 
by the subinterval, (u — K 1 , u + Ki), while to ensure Tu is onto Q* we need only assume 
that f2 is G-convex with respect to {y G il*,z G I such that |G(-,y, z) — u | < Kx}. 
Furthermore if (4.6) holds then we can take any uq > mo in Corollary 4.6. 

Theorems 4.2,4.5 and Corollary 4.6 clearly apply to the parallel beam example (1.18) 
and under the hypotheses of Corollary 4.6 we would obtain the existence of a positive 
elliptic solution u of equation (1.18) with \Du\ < 1. 

We conclude this section with some remarks about the point source reflection problem 
treated in [10, 11, 17]. Invoking the ellipsoid of revolution in its polar form [10, 11] 
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and projecting onto the unit ball in W 1 , we may take U = -Bi(O) x R n , / = (0, oo). In 
particular, if the target hypersurface lies in a graph y n+1 = T(y) in R™ +1 , we may model 
the reflection process by a generating function: 

(4.8) G(x, y, z) = l -y z +\y\> + Ti- x.y - y/l - \x\* t}. 

Differentiating we obtain, 

(4-9) Gx = -(-y+ i , T), 

z Jl — \x\ 2 

G ^ = —R ) iono/o (/ + x <S> x)T. 

z(l — \x\ I ) 6 l 2 

Let us restrict here to a special case where the target lies in a hyperplane not above the 
source at the origin, namely y n =\ = r < . Writing u = u — p.x,p = G x , it follows that 
u > \p\ and moreover we have the formulae: 

(4.10) Z=(l- , 2 ™ )/(** - \p\% 



Y = -pZ + 



\x\ 

TX 



X 



2 



A = v . = (I + x®x). 

(1 - \x\ 2 )(Jl- \x\ 2 -2tu) 



It is now easy to check that G satisfies G1,G2,G1*,G2,G3,G4 with I(x,y) = J(x,y) = 
(0, oo) if r < 0, and Gl,G2,Gl*,G3w,G4w for r = 0, in which case A = 0. Furthermore 
if Q G -Bi(O), then G also satisfies (4.6). The resultant existence and regularity results 
correspond to special cases of [11] except we are also able to solve the boundary value 
problem (1.7) in the classical sense when Q is also G-convex. However the local regularity 
results in [11] are more general in that condition G 3 need not be satisfied everywhere as 
in this case a stronger version of Lemma 2.2 enables local regularity to be proved in the 
sub-domains where G3 is satisfied. When r = 0, G-convexity coincides with convexity 
and local regularity under G*-convexity of the target domain fl* follows from [1] , Lemma 
3. 

Finally we note that the existence of globally smooth solutions of the point source 
and parallel beam near field reflection problems is established in [17] and that solutions of 
these and related problems can also be represented as potentials of a nonlinear Kantorovich 
problem, [7, 15, 16]. 
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